Subsurface Structure Research Articles

Dermoscopic structures and patterns used in melanoma detection.

Melanoma is the leading cause of skin cancer-related deaths. Yet, early detection remains the most cost-effective means of preventing death from melanoma. Early detection can be achieved by a physician and/or the patient (also known as a self-skin exam). Skin exams performed by physicians are further enhanced using dermoscopy. Dermoscopy is a non-invasive technique that allows for the visualization of subsurface structures that are otherwise not visible to the naked eye. Evidence demonstrates that dermoscopy improves the diagnostic accuracy for skin cancer, including melanoma; it decreases the number of unnecessary skin biopsies of benign lesions and improves the benign-to-malignant biopsy ratio. Yet, these improvements are contingent on acquiring dermoscopy training. Dermoscopy is used by clinicians who evaluate skin lesions and perform skin cancer screenings. In general, under dermoscopy nevi tend to appear as organized lesions, with one or two structures and colors, and no melanoma-specific structures. In contrast, melanomas tend to manifest a disorganized pattern, with more than two colors and, usually, at least one melanoma-specific structure. This review is intended to familiarize the reader with the dermoscopic structures and patterns used in melanoma detection.

Inclusion of bedrock vadose zone in dynamic global vegetation models is key for simulating vegetation structure and function

Abstract. Across many upland environments, soils are thin and plant roots extend into fractured and weathered bedrock where moisture and nutrients can be obtained. Root water extraction from unsaturated weathered bedrock is widespread and, in many environments, can explain gradients in vegetation community composition, transpiration, and plant sensitivity to climate. Despite increasing recognition of its importance, the “rock moisture” reservoir is rarely incorporated into vegetation and Earth system models. Here, we address this weakness in a widely used dynamic global vegetation model (DGVM; LPJ-GUESS). First, we use a water flux-tracking deficit approach to more accurately parameterize plant-accessible water storage capacity across the contiguous United States, which critically includes the water in bedrock below depths typically prescribed by soil databases. Secondly, we exploit field-based knowledge of contrasting plant-available water storage capacity in weathered bedrock across two bedrock types in the Northern California Coast Ranges as a detailed case study. For the case study in Northern California, climate and soil water storage capacity are similar at the two study areas, but the site with thick weathered bedrock and ample rock moisture supports a temperate mixed broadleaf–needleleaf evergreen forest, whereas the site with thin weathered bedrock and limited rock moisture supports an oak savanna. The distinct biomes, seasonality and magnitude of transpiration and primary productivity, and baseflow magnitudes only emerge from the DGVM when a new and simple subsurface storage structure and hydrology scheme is parameterized with storage capacities extending beyond the soil into the bedrock. Across the contiguous United States, the updated hydrology and subsurface storage improve annual evapotranspiration estimates as compared to satellite-derived products, particularly in seasonally dry regions. Specifically, the updated hydrology and subsurface storage allow for enhanced evapotranspiration through the dry season that better matches actual evapotranspiration patterns. While we made changes to both the subsurface water storage capacity and the hydrology, the most important impacts on model performance derive from changes to the subsurface water storage capacity. Our findings highlight the importance of rock moisture in explaining and predicting vegetation structure and function, particularly in seasonally dry climates. These findings motivate efforts to better incorporate the rock moisture reservoir into vegetation, climate, and landscape evolution models.

Solutions of Surface-Wave Dispersion and Attenuation in Stratified Viscoelastic Media Using a Spectral-Element Approach

ABSTRACT Efficient and accurate calculation for the dispersion and attenuation of the surface waves in viscoelastic media is numerically challenging because the eigen wavenumbers are located in the complex domain. In this study, we propose a semianalytical spectral-element method (SASEM), which can determine the complex eigen wavenumbers by solving linear eigenvalue problems. By simplifying the structure of the eigenvalue problem, we significantly improve the calculation efficiency. The implementation of the frequency-dependent automatic discretization, semi-infinite element, and mode filter guarantees the correctness and accuracy of the modal solutions. Because no root-searching schemes are required, the root-skipping problem is naturally avoided. The numerical tests show that the SASEM can provide sufficiently accurate solutions with much less computation cost than traditional Muller’s method. Meanwhile, SASEM exhibits high flexibility when applied to media the parameters for which vary continuously with depth. To demonstrate the effectiveness of SASEM for complicated dispersion features, the dispersion curves and eigen wavefields of the viscoelastic media with a low-velocity layer are also analyzed. The results of numerical tests indicate the versatility, efficiency, and accuracy of our method. With further study, the proposed SASEM has the potential to become a promising tool for the investigation and retrieval of viscoelastic subsurface structures.

Temperature Dependence of the Dielectric Constant on the Lunar Surface Based on Mini-RF and Diviner Observations

Radar observation is an effective way to understand subsurface structures in terms of the dielectric constant, whose controlling factors include chemical composition, packing density, and water/ice content. Recently, laboratory measurements have shown that the dielectric constant of lunar regolith simulants also depends on the temperature, which has never been evaluated from remote sensing data. In this study, we estimated the dielectric constant from the Miniature Radio Frequency (Mini-RF) data on a lunar crater floor in the north polar region at two different local times (i.e., different surface temperatures). We calculated the dielectric constant using the inversion method and obtained the bolometric surface temperature from the Diviner Lunar Radiometer Experiment (Diviner) data. The histograms of the estimated dielectric constant values are different between the two local times. This could be interpreted as a result of the temperature dependence of the dielectric constant, while further evaluation of the influence of topography on the incidence angle and small surface roughness is needed. Nevertheless, our result suggests that the temperature dependence of the dielectric constant should be considered when interpreting S-band radar observations of the Moon and other celestial bodies with large surface temperature differences.

A depth-variant seismic wavelet extraction method for basis pursuit inversion with an impedance trend constraint

The seismic images produced by prestack depth migration indicate more accurate subsurface structures than time images, resulting in a growing need for depth-domain inversion. However, due to the strong nonstationarity exhibited by depth-domain seismic data, time-domain inversion methods based on the convolutional model cannot be directly applied in the depth domain. To address this issue, we develop a method for extracting a depth-variant seismic wavelet, which is then combined with a nonstationary convolutional model to enable direct inversion of the depth-domain acoustic impedance (AI). First, we extend the Morlet wavelet to the depth domain and develop an orthogonal matching pursuit spectral decomposition method using the depth-domain Morlet wavelet. We then investigate the waveforms and wavenumber spectra similarities between the depth-domain Morlet wavelet and depth-domain Ricker wavelet and extract depth-variant Ricker wavelets from the depth-wavenumber spectrum. We add a depth-domain impedance trend constraint to the conventional basis pursuit inversion to enhance the lateral continuity of the inversion results. Then, we attain direct inversion of the depth-domain AI. Tests of synthetic and field data demonstrate that our method achieves high-accuracy inversion results while maintaining high computational efficiency, highlighting our approach’s effectiveness and strong reservoir characterization potential.

JOINT INVERSION OF DIRECT WAVES, PRIMARIES AND MULTIPLES FOR VERTICAL CABLE SEISMIC DATA

By anchoring the vertical cable in seawater and receiving the data excited by the towing cable seismic source, the reflection data from subsurface structures with a high signal-to-noise ratio, wide frequency band, and wide azimuth angle can be obtained, thereby achieving high-quality imaging of the target area. However, due to drifting receivers driven by ocean currents and nonuniform coverage caused by the vertical orientation of the receivers, the traditional imaging method cannot be directly used in vertical cable seismic (VCS) data. Furthermore, crosstalk between primaries and multiples commonly contaminates the final imaging. A new inversion method, based on the Bayesian inference and the assumption that all errors are Gaussian distributed, is presented to address the above problems. The direct waves are used to invert source/receiver locations. The primary and multiple wavefields are accurately reconstructed based on the reconstructed direct wavefield, and the corresponding reflectivity updates are obtained. The model-domain Hessian matrices consisting of different data are approximated by their diagonal elements to mitigate the nonuniform coverage issue, and the related covariance matrices and data-domain Hessian matrices are defined and approximated, respectively, to alleviate the calculation problem. In addition, the zero-lag correlated wavefields in the gradient are consistent with the actual wavefield propagation process, e.g., a primary wave is generated when the direct wavefield encounters reflectivity, avoiding the crosstalk. Numerical examples for synthetic models and field data illustrate that the proposed method can accurately recover the source/receiver locations and provide high-resolution, balanced-amplitude imaging results without evident crosstalk.

Introducing a tailored site delineation approach to optimize the design of managed aquifer recharge surface spreading infrastructure

Managed aquifer recharge (MAR) is an emerging solution to effectively replenish overused groundwater resources, but high associated costs often hinder its uptake. There are several parameters to consider when determining the cost of MAR, including recharge types, MAR schemes scale, land acquisition, operation and maintenance period, and the hydrogeological setting. Hydraulic conductivity and its spatial variability are the most significant hydrogeological parameters for predicting infiltration rates at MAR sites, but limited data availability makes accurate predictions often challenging. Hence, it is essential to increase data availability to optimize MAR efficiency and to better assess its economic performance. Therefore, a novel approach for MAR planning is presented in this article to efficiently enhance the data availability during the conceptual stage of MAR planning by combining site exploration, spatially resolved infiltration rate estimation, and cost analysis. The approach was applied at an actual MAR site in Vincenza, Italy, and incorporated an electromagnetic induction survey and direct push profiling as across-scale investigation methods to enhance site delineation and subsequent parameterization of the identified zones using laboratory analysis of collected samples. This MAR site is particularly suitable for its heterogeneous subsurface structure and potential to counteract groundwater depletion. The performed investigation resulted in an image of the lateral variation of subsurface conditions and allowed the delineation of three zones with different infiltration behavior, with 99.5% of infiltration occurring in just two zones, representing only 22.1% of the site. In this particular scenario, the overall cost per cubic meter of recharged water can be reduced by 59.1% by identifying and eliminating unfavorable zones. This study provides scientists and practitioners with a useful tool for MAR planning that can be applied to a wide range of sites with complex geology, thereby reducing water costs, minimizing the environmental impact of infrastructure construction, and reducing land use conflicts.

Stochastic seismic acoustic impedance inversion via a Markov-chain Monte Carlo method using a single GPU card

Seismic acoustic impedance inversion plays an important role in understanding subsurface structures and obtaining subsurface properties. The stochastic approach is one of the methods used for impedance inversion, and it aims to produce more reliable results by accounting for modeling uncertainty. Stochastic inversion represents the uncertainty of a subsurface model as a probability distribution and uses this distribution to estimate model parameters. In this study, seismic acoustic impedance inversion was performed using a Markov-chain Monte Carlo approach, which is a sampling method used in the stochastic process. Utilizing Bayesian inference based on prior information and observed data, we implemented acceptance probabilities and performed acoustic impedance inversion for field data through iterative calculations. We employed a single GPU card to execute the inversion algorithm and conducted a comparative analysis with the results obtained using a cluster composed of multiple CPU cores. Through a computational speed analysis, the efficiency of the algorithm using a GPU was verified, while uncertainty analysis was employed for algorithm validation. Through these analyses, we confirmed the feasibility of applying our developed GPU algorithm to tasks that require the inversion of extensive data, such as high-resolution seismic exploration and CO2 storage monitoring.

Determination of Hydrocarbon Production Zones using Well Logging Data in the "North_Rawa" Well of the South Sumatra Basin

Petrophysical analysis with well logging is carried out to determine the value of petrophysical parameters of rocks in the formation. Then to find out information about the subsurface structure as well as the depth zone of layers containing hydrocarbons and the types of hydrocarbons from the production field. In this research, an evaluation of reservoir formation was carried out in the North_Rawa field of the South Sumatra Basin, based on petrophysical log analysis from data from three wells in the North_Rawa field. Reservoir identification results from physical rock parameters such as porosity, shale content, water saturation and permeability obtained from formation evaluation calculations. Porosity is calculated to determine the pores, shale content is calculated to determine the volume of shale in the reservoir, water saturation is calculated to determine the water content in the formation, and permeability is calculated to determine the rock's ability to pass fluid. The results of petrophysical calculations were obtained in the hydrocarbon prospect zone in the North- rawa 1 well at a depth of 3387.5 – 3938 feet, the North- rawa 3 well at a depth of 2667.5 – 3735.5 feet, the North- rawa 4 well at a depth of 2759.5 – 3775 feet. The average porosity values for each well are 8%, 15.4%, 16%. The average shale content value for each well is 61%, 48%, 53%. Average water saturation values are 61%, 54%, 68%. Average permeability values 46.7 mD , 98.9 mD , 3.38 MD. The hydrocarbon prospect in the research well contains natural gas and oil hydrocarbons with sandstone and shale formation lithology.

Improved Approaches for 3D Gravity and Gradient Imaging Based on Potential Field Separation: Application to the Magma Chamber in Wudalianchi Volcanic Field, Northeastern China

The gravity and gradient anomalies contain valuable information about the underground geological structures at various depths. Deep and shallow buried source bodies are able to be identified through multi-scale field separation processes, and visual comprehensions of geological structures can be obtained via 3D density inversion techniques. In this study, we propose an improved 3D imaging strategy based on gravitational field separation using the preferential continuation filter. This strategy incorporates the relationship between spectral features and buried depths of source bodies, allowing for a one-step transformation from planar gravity and full-tensor gradient field observations to a 3D density structure in the wave-number domain. Synthetic tests validate the effectiveness and robustness of the gravity and gradient imaging approaches, highlighting their advantages in high vertical resolution and low computational requirements. Nonetheless, it should be noted that the imaging effects of horizontal gradients Γxx and Γyy are unsatisfactory due to their weak noise resistance. Thus, they are not suitable for real data applications. The other imaging approaches are further applied to recover the subsurface 3D density structure beneath the Weishan cone in Wudalianchi Volcanic Field, Northeastern China. Our results provide insights into the possible location and shape of the low-density magma chamber. Also, the potential presence of partial melts is inferred and supported from a gravity perspective. The primary advantage of these approaches is their ability to generate a reasonable geological model in scenarios with limited prior information and physical property constraints. As a result, they have significant practical value in the field of applied geophysics, including mineral exploration and volcanology studies.

Reliable Tilt-depth estimates based on the stable computation of the tilt angle using robust vertical derivatives

The Tilt-depth is a popular approach for determining depths of magnetic sources. As this method is based on the distance between contour levels of the tilt angle, it can lead to inaccurate depth estimates when the contour lines are distorted by the presence of noise. In this paper, we stabilize the Tilt-depth method based on the computation of stable vertical derivative obtained by the β-VDR method. The presented method is demonstrated on synthetic magnetic anomalies and real magnetic data from the Arabian Shield. The results obtained from the synthetic examples coincide well with the actual depths. These results proved the utility of the presented method in cases where the field is corrupted by noise. The real example shows that the presented method can provide valuable information on subsurface structures of the area where the Tilt-depth estimates are consistent with the result of the horizontal tilt angle. The findings show that the presented method is less sensitive to noise and can provide source edges and depths more clearly and with higher accuracy.

Application of Generalized Reciprocal Method on 2D Seismic Refraction Data in Mt. Manglayang, West Java

The imaging of shallow subsurface structures, weathered rock thickness and velocity propagation distribution of the rocks can be identified by seismic refraction. This method is one of the geophysical exploration methods utilizing refracted wave once it reaches the boundary of subsurface layer. In this research we used the generalized reciprocal method (GRM) as one of the robust processing methods in analyzing subsurface data. This method was chosen due to its accuracy in interpreting shallow subsurface layer with highly undulating refractors by determining time velocity analysis, XY optimum distances, and time depth analysis, then, the expected depth values can be achieved. The acquisition of data for this research was conducted using 13 geophones with forward and reverse sources, the data were then picked to get travel time values and inverted to obtain real geological setting of the earth. The results were interpreted as 2 layers, the first layer had a velocity distribution of 499.289 m/s which was identified as a weathered layer with a thickness of about 4.74 meters, whereas the second layer was interpreted as clay rock with velocity distribution of 1270.433 m/s with the thickness reached up to 16 meters.
 Keywords: generalized reciprocal method, 2D seismic refraction

Subsurface architecture of the Wajrakarur Kimberlite Field, Dharwar Craton, India: Evidence from ground magnetic surveys

The Ground magnetic study in the Wajrakarur Kimberlite Field and adjoining area has unveiled prominent geological variations, including felsic doming beneath the PGC-II formation to the west of the Cuddapah Basin and mafic doming beneath the Closepet Granite zone. These doming features within the upper mantle are identified as sources of mantle upwelling, which influence the emplacement of corresponding intrusive rocks. The presence of mafic chambers enhances the magnetic response, leading to localized high signatures in magnetic data and contributes to a better understanding of the subsurface composition and distribution of rock types. The present magnetic data throw lights to study the complexity of the sub-surface structure and geological activities, have shaped the study area and serve as valuable markers for understanding the geological evolution of the area. The intersections of magnetic lineament oriented in NW-SE, NE-SW, E-W yield crucial insights into geological structures and the emplacement of Kimberlite intrusions. Through comprehensive magnetic data analysis, the study offers a unique cross-sectional view of the subsurface, reaching depth up to 60 km. This perspective unveils magnetic sources and structures, revealing shallow features, interfaces of the Archaean basement, crustal roots, and upper mantle configuration. Profound magnetic gradients observed in the study area serve as indicators of fault lines influenced by tectonic activity, kimberlite intrusion, and other fault-related mineralization if any. These features contribute to a deeper understanding of the structural configurations and geological history of the area. The present study provides information about the lateral and vertical layouts of major structures as well as the deep-seated channels. It also introduces the mafic and felsic categories of rocks, explaining the complex relationships between underlying structures, geological characteristics, and magnetic anomalies.

Characterization of subsurface geology and hydrogeology in Kribi -Cameroon using electrical resistivity soundings and 3D-Implicit modelling: Baseline for groundwater resource management

This scientific article aims to characterize the geology and hydrogeology of Kribi, Cameroon for the purposes of water resources management, infrastructure development, and environmental protection. The specific objectives include mapping the subsurface electrical resistivity structure, identifying groundwater potential zones, and constructing a hydrogeological model for groundwater modelling. The study employed the Schlumberger method to map the electrical resistivity structure of the 3D subsurface and utilized 3D implicit modelling techniques to delineate groundwater potential zones. The study revealed that Kribi exhibits a highly complex geology, characterized by three-layered electrical resistivity models, typical of sedimentary and metamorphic formations. The major lithologies identified include laterite, clay, sand, quartz, and granite. The region encompasses areas with low water potential and brackish zones alongside high groundwater production areas. The study established the presence of aquifers within Kribi, with the unconfined aquifer found at depths ranging from 2 m to 13 m, and the confined aquifer extending from 14 m to 40 m. The aquifers are notably thick and productive, indicating substantial groundwater resources. Static groundwater volumes were estimated at 50,179,381 m3 and 317,118,946 m3 for the unconfined and confined aquifers respectively. Despite this potential, water resources management in the area is non-existent. The study recommends the implementation of Integrated Water Resources Management (IWRM) to improve water resource management, with a primary focus on groundwater modelling. The northern and southern regions of Kribi exhibit high permeability, which is advantageous for groundwater recharge and storage projects but also poses challenges for environmental remediation efforts. Consequently, safeguarding the delineated area for groundwater provision is of utmost importance.

Near surface resistivity investigation of regolith of Middle Benue Trough: implications on road pavement along Lafia-Shendam Axis, North Central Nigeria

ABSTRACT Electrical characteristics and subsurface structures associated with the regolith of Middle Benue Trough have been investigated. Vertical electrical sounding (VES) and constant separation traversing (CST) techniques were employed for data acquisition along Lafia-Shendam Road, Lafia, North Central Nigeria. Forty (40) Schlumberger sounding stations and four (4) Wenner resistivity profiles were occupied along opposite sides of the road to delineate and decipher the subsurface geological structures. Second-order Dar Zarrouk parameters (longitudinal and transverse resistivities, mean resistivity and coefficient of anisotropy) computed from VES parameters were also used to evaluate the pseudo-anisotropy characteristics of the regolith cover. Geoelectric sequences delineated within the study area include thin topsoil, alluvium deposit, weathered layer (laterite and saprolite) and fresh basement. The regolith is hosted by a weathered sedimentary rock at an average depth which varies from 10.2 to 37.3 m. The regolith structure is not complex and has a greater proportion of conductive laterite. The saprolite associated with the host sedimentary rock, masked in the VES technique, is revealed as a very thin zone with moderate resistivity in the CST 2-D structures. The predominance of laterite in the regolith is also corroborated by the range of values of the coefficient of anisotropy and mean resistivity.

Integrated geophysical methods to constrain subsurface structures of Tulu Moye-Bora-Berecha axial volcanic complex, main Ethiopia rift: Implications for geothermal resources

The Main Ethiopian Rift (MER) is a well-known continental rift whose axial sector is characterized by the occurrence of regularly spaced silicic caldera complexes and central stratovolcanoes, interspersed with large fields of fissural basalts, small mafic scoria cones and numerous young normal faults and fissures. The Tulu Moye-Bora-Berecha volcanic complex is found in the central portion of the MER and includes the Tulu Moye geothermal prospect area. A combination of gravity and magnetic methods was used to better constrain the subsurface volcanic stratigraphy and tectonic structures. Regional and residual anomaly maps were produced from the gravity data and the magnetic data were corrected to produce anomaly and enhanced maps. A complete Bouguer anomaly contour map was produced after the necessary reduction was applied to the gravity data. Due to the geomagnetic field's dipolar nature and since the study area is in the equatorial region (<15°), the reduced-to-equator (RTE) technique was used to minimize external effects and correct the data as if the body had been laid at the magnetic equator. The anomalous source's depth was calculated using a Euler depth solution and spectral analysis approach. Joint 2D forward models on three profiles were developed using GM-SYS of Oasis Montaj software. From the interpretation of the geophysical results, the following conclusions have been reached: (1) the crystalline basement is more raised around Salen ridge and west of it; (2) the main heat source of the geothermal system appears to be a central region of the studied area near Salen ridge and is estimated to be at 4–5 km depth, (3) under Gnaro obsidian dome, the basaltic and silicic volcanic horizon is thin, and the deep-seated regional fault serves as the main conduit for the passage of hot fluid to the surface and (4) the gravity and magnetic anomaly plots, regional-residual maps and enhanced data plots all indicate towards the existence of a major geological feature-a large caldera-comprising of Tulu Moye, Bora and Berecha volcanic centers.

Subsurface structure of Bali Island inferred from magnetic and gravity modeling: new insights into volcanic activity and migration of volcanic centers

Subsurface structure of Bali Island inferred from magnetic and gravity modeling: new insights into volcanic activity and migration of volcanic centers

Improving full-waveform inversion based on sparse regularisation for geophysical data

Abstract Full waveform inversion (FWI) is an advanced geophysical inversion technique. FWI provides images of subsurface structures with higher resolution in fields such as oil exploration and geology. The conventional algorithm minimises the misfit error by calculating the least squares of the wavefield solutions between observed data and simulated data, followed by gradient direction and model update increment. Since the gradient is calculated by forward and backward wavefields, the high-accuracy model update relies on accurate forward and backward wavefield modelling. However, the quality of wavefield solutions obtained in practical situations could be poor and does not meet the requirements of high-resolution FWI. Specifically, the low-frequency wavefield is easily affected by noise and downsampling, which influences data quality, while the high-frequency wavefield is susceptible to spatial aliasing effects that produce imaging artefacts. Therefore, we propose using an algorithm called sparse relaxation regularised regression (SR3) to optimise the wavefield solution in frequency domain FWI, which is the forward and backward wavefield obtained from the Helmholtz equation, thus improving the FWI's accuracy. The sparse relaxation regularised regression algorithm combines sparsity and regularisation, allowing the broadband FWI to reduce the effects of noise and outliers, which can provide data supplementation in the low-frequency band and anti-aliasing in the high-frequency band. Our numerical examples demonstrate the wavefield optimisation effect of the sparse relaxation regularised regression-based algorithm in various cases. The improved algorithm's accuracy and stability are verified compared to the Tikhonov regularisation algorithm.

Automatic Extraction of Martian Subsurface Layer from Radargrams Based on PDE Denoising and KL Filter

The polar regions of Mars, including the South and North Poles, are crucial for studying Martian climate and geological history, as they contain the largest reservoir of subsurface water ice. This study introduces a new approach for reflector detection, which includes radargram denoising to effectively enhance the signal of underground reflectors, peak detection to extract the positions of subsurface stratification from the radar echoes, and peak points connection to form continuous layers. The mapped enhancement denoising process involves a linear brightness adjustment and a fourth-order diffusion equation to enhance the signal of the subsurface layers for effective detection. The subsurface detection extracts the surface and subsurface peak points based on a peak detection algorithm, while using locally window-enhanced peak filtering and Kullback–Leibler (KL) divergence mapping to filter out non-stratified peak points. Finally, the layered connection process uses the proximity parameter to connect peak points in the same layer. Applied to multiple SHARAD (Shallow Radar) images at the Martian poles, this algorithm demonstrated a false detection rate below 5%. Compared to other methods, this method has a missed detection rate of less than 5% and, additionally, exhibits fewer discontinuities in layer connectivity. Therefore, this algorithm shows exceptional proficiency and applicability in analyzing the complex subsurface structures of the Martian polar regions.

Integrated stratigraphical, hydrogeochemical, and geophysical studies to characterize the shallow groundwater aquifers in Wadi Araba area, western Gulf of Suez, Egypt

The construction of the New Galala City besides great mining (e.g., konoozmisr) and Seawater Reverse Osmosis Desalination Plant projects in the semiarid Gulf of Suez region of Egypt increases the demand for groundwater resources. So, the current research provides integrated stratigraphical, hydrogeochemical, and geophysical studies to investigate the shallow groundwater aquifers in Wadi Araba area on the western side of the Gulf of Suez. The stratigraphic column starts with pre-Cambrian basement rocks which are non-conformably overlain by sedimentary succession of Phanerozoic Erathem. A detailed land magnetic profile has been conducted and processed to estimate the aquifer’s geometrical parameters. This includes the subsurface spatial extension and thickness to delineate the prevailed subsurface structures in the study area. Lithostratigraphically, these aquifers are represented by the Paleozoic- lower Cretaceous siliciclastic-dominated sediments of the Nubia Group and fractured carbonate-dominated facies of the upper Cretaceous Galala Formation. Lithological characters, sedimentary facies (fluvial to shallow and open marine), microfacies types, and diagenesis processes (e.g., compaction, cementation, dissolution, and fracturing) are the main factors that affect the quality of the shallow aquifers in the study area. The main microfacies of the studied aquifers are quartz arenite, subarkose, sublitharenite, mudstone/wackestone, and packstone/grainstone. These facies are characterized by different primary and secondary porosities. According to the water sampling standard methods, nineteen groundwater samples were collected from the available hand-dug wells and springs from the Nubian sandstone aquifer and the Upper Cretaceous aquifer, respectively. Inductively coupled plasma–mass spectrometry (ICP–MS), titration, spectrophotometer and a flame photometer were used to obtain a complete chemical analysis. The results of hydrochemical analyses indicated that total dissolved solids (TDS) and concentrations of major ions as well in the lower Cretaceous aquifer increase toward the Gulf of Suez and in general higher than that of springs of the Upper Cretaceous Aquifer. The predominant factors influencing groundwater chemistry include seawater intrusion, evaporation, mineral alteration, dissolution of fossils, and ion exchange. Groundwater in the study area is unsuitable for drinking except in one sample and moderately suitable for irrigation with caution for higher salinity in samples located near the Gulf of Suez.